Photon recycling in Fabry-Perot micro-cavities based on SisN 4 waveguides. 
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We present a numerical analysis and preliminary experimental results on one-dimensional Fabry- 
Perot micro-cavities in SiaN4 waveguides. The Fabry-Perot micro-cavities are formed by two dis- 
tributed Bragg reflectors separated by a straight portion of waveguide. The Bragg reflectors are 
composed by a few air slits produced within the SiaN4 waveguides. In order to increase the quality 
factor of the micro-cavities, we have minimized, with a multiparametric optimization tool, the in- 
sertion loss of the reflectors by varying the length of their first periods (those facing the cavity). To 
explain the simulation results the coupling of the fundamental waveguide mode with radiative modes 
in the Fabry-Perot micro-cavities is needed. This effect is described as a recycling of radiative modes 
in the waveguide. To support the modelling, preliminary experimental results of micro-cavities in 
SisN4 waveguides realized with Focused Ion Beam technique are reported. 
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I. INTRODUCTION 



Silicon-based photonics is the key technology for 
manipulating, controlling, and detecting light at sub- 
micrometer length scales in silicon 0. Owing to the 
high index contrast between silicon based materials and 
air, these systems are ideal to study devices associated 
to photonic band-gap materials. Among these, electro- 
magnetic resonant cavities which are able to trap light 
can be considered as building blocks of future photonic 
circuits 0. In these systems, the confinement of the 
photons within a finite volume is assured by a periodic 
refractive index modulation of the surrounding medium. 
The best choice is a three-dimensional refractive index 
modulation that assures a full confinement of photons 
in the cavities. Nevertheless, there are lot of difficulties 
in fabricating three-dimensional periodic structures op- 
erating at infrared wavelengths. As a consequence, it 
seems to be favorable to explore new devices with two 
or one dimensional (ID) refractive index modulation for 
which the fabrication technology is well established; light 
confinement in the other dimensions can be realized by 
total internal reflection within optical waveguides. Pho- 
tonic crystals based on such an approach exhibit a quasi- 
photonic band-gap due to the lack of three dimensional 
confinement. An implementation of ID photonic crystals 
are photonic crystal slab waveguides, where a high-index 
core layer is sandwiched between lower index claddings. 
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The ID refractive index modulation is achieved by pro- 
ducing air trenches across the channel waveguide. A 
ID photonic crystal with a defect is also described as 
a Fabry-Perot micro-cavity. This is formed by two first 
order Bragg mirrors separated by a spacer. The Bragg 
mirrors are constituted by a sequence of A/2 thick peri- 
ods, where each period contains an air slit and a waveg- 
uide segment. A photon mode propagating in the channel 
waveguide excites many cavity radiative modes when it 
is transmitted trough the cavity. This fact degrades the 
quality factor (Q-factor) of the cavity. Indeed the key 
point for high Q-factor cavities relies on a fine tuning 
of the first periods facing the cavity spacer 0, 0, 13 to 
reduce the excitation of radiative cavity modes and de- 
crease the impedance mismatch between waveguide and 
cavity modes. Nevertheless the physics beyond these ef- 
fects is not completely understood, and it is still under 
debate 

The aim of this work is to characterize and optimize the 
performance of a ID Fabry-Perot (FP) micro-cavity cen- 
tered at 1550 nm and formed in a SisN 4 channel waveg- 
uide. We start our characterization optimizing the in- 
sertion losses of a Distribute Bragg Reflector (DBR) via 
engineering the first periods facing the waveguide spacer, 
as proposed in 0. Then we use the optimized mirrors 
to build a FP micro-cavity and we study the quality fac- 
tor (Q-factor) of the optimized system as a function of 
the resonance order. In the second part of the work we 
show that the classical FP model of Q-factor is unable 
to explain the numerical results. These can be under- 
stood by considering the recycling of radiative modes 
in the waveguide 0. In the last part of this work, we 
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Figure 1: a) Lateral view of the simulated system composed 
by the input waveguide and the mirror where the coupling 
periods are emphasized, b) Insertion losses of six optimized 
mirrors, L=l-R-T, as a function of the wavelength. The mir- 
rors differ by the number of coupling periods which face the 
input waveguide. The arrow is a guide for the eye to follow 
the results as the number of coupling periods increases. 

report preliminary experimental data on a FP micro- 
cavity. The FP micro-cavity has been realized with Fo- 
cused Ion Beam (FIB) processing, starting from a single 
mode S13N4 channel waveguide and removing waveguide 
slices to define the photonic structure. Our modelling is 
able to explain the experimental data. The full theoret- 
ical characterization and realization of other devices are 
currently in progress to further validate the numerical 
predictions. 

II. ENGINEERED MIRROR FOR HIGH 
Q-FACTOR DEVICE 

Classical FP micro-cavities on waveguides are formed 
by a straight waveguide surrounded by two identical 
DBR. The main limitation of these systems is the modal 
mismatch between the fundamental cavity mode (waveg- 
uide mode) and the fundamental Bloch mode of the DBR. 
This mismatch causes strong coupling of the resonant 
guided mode with radiative modes of the DBR, decreas- 
ing the performance of the micro-cavity @ . Thus an op- 
timization of the coupling between waveguide and DBR 
is needed. The first step is to consider the problem re- 
lated to the reflection of a transverse electric (TE) mode 
impinging from a monomodal waveguide onto a DRB (see 
Figure ||TJi) for a sketch). 

As material system we choose S13N4 waveguides since 



we are investigating them as substitute of more expen- 
sive silicon on insulator waveguide 0- A monomode 
SiaN4 waveguide formed by a multilayer core of SiaN4 
and SiC>2 and a top-bottom cladding of Si02, will be 
considered hereafter (see the experimental section Hvl for 
more details) . The waveguide is interfaced to a first order 
Bragg mirror with a period of 490 nm which is composed 
by air slits (width 100 nm) etched 1/im down through 
the bottom cladding. With these parameters the DBR 
supports a single guided TE polarized Bloch mode - a 
mode below the light-line. The characterization of the 
waveguide-mirror interface is made by exciting the fun- 
damental waveguide mode and calculating the reflection 
R (transmission T) coefficient of this mode through the 
waveguidc-DBR system. The insertion loss spectra, de- 
fined as L = 1 — R — T, are calculated from 1.4 fim 
to 1.7 /im, which includes the DBR stop band region. 
All the simulations have been made with a commercial 
software flOl based on an eigenmode expansion (EME) 
method |ll|. The minimization of the insertion losses is 
made by means of a multiparametric optimization tool. 

Optimization of the insertion losses is done by inserting 
coupling periods at the interface between the waveguide 
and the DBR. The physical lengths of each air slit and 
waveguide segment in the coupling periods are the free 
parameters of the multiparametric optimization. The nu- 
merical results of the optimization are shown in Figure 
||TJ)) where the insertion loss spectra of six optimizations 
are reported as a function of wavelength. The six curves 
correspond to different numbers of coupling periods, var- 
ied from zero (no coupling periods) to five. With in- 
creasing the number of coupling segments a well defined 
minimum appears at 1.55 /im; the insertion losses are 
decreased by two order of magnitude, from 10 -2 for the 
bare mirror to 3 x 10 -4 for the five periods engineered 
mirror. Table 1 reports the parameters of the coupling 
periods for the five mirrors. For each engineered mirror 
we report the period length A and the filling fraction f.f. 
(the ratio between the length of air slit and the period) 
of each set of coupling periods. 

The engineered mirrors have been then used to form 
optimized micro-cavities (Figure [5J|. Figure [21 reports 
their calculated Q-factor as a function of the resonance 
order m when no (circles), one (triangles) and two 
(squares) coupling periods are used. The micro-cavity 
intrinsic Q-factor was calculated from the asymptotic 
value of the full width at half maximum AA of the reso- 
nant peak as the number of DBR periods increases while 
the number of coupling periods is kept constant; know- 
ing the micro-cavity resonance wavelength Ao = 1.55/iin, 
Qint = Ao/AA. The resonance order is defined as in 
Eq. |2| For the smallest defect length the minimum value 
of m is about 4, due the large penetration depth of the 
mode defect in the mirror (see also the discussion after 
Eq.C|. 

It is seen that the Q-factor increases as the number 
of coupling periods increases, which was expected, since 
the insertion losses of the mirrors L decrease as shown 
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Tablel: period length A and filling fraction f.f. of each set 
of coupling segments for the five engineered mirror. For the 
periodic mirror A = 490nm and f.f.=0.204. 



in Figure QJ. For instance at m rj 9 we have Q w 3 x 
10 3 , 1.1 x 10 4 , 5.1 x 10 4 for no, one and two coupling 
periods, respectively. On the other hand, it is clearly seen 
that such an increase is not very regular as m changes. 
Such a behavior is discussed in the next section. 



III. INTERPRETATION OF THE SIMULATION 
RESULTS 

Let us consider a FP cavity, with mirrors characterized 
by complex reflectivity and transmissivity coefficients, 
r = |r|e !¥,m and t = \t\e lVm , respectively; the spacer 
between the two mirrors is characterized by a physical 
length L and by an effective refractive index n e ft. The 
expression of the cavity transmission is given by the Airy 
formula 



tFP 



tV* 



1 



(1) 



where the phase accumulated by the electromagnetic 
wave after half round-trip in the cavity is 



* = —n e ffL + tp„ 



(2) 



Under the assumption that 1 - 



<C 1, the quality factor 



of the micro-cavity can be expressed as 



Q = 



l-|rf 



27T dlfm 



x=x„ 



1 - |r| : 



r m7r (3) 



where Ao is the resonant wavelength, n g — n e ff — 
\(dn e ff/d\) is the group index of cavity mode. In the 
eq. © we have also identified the resonance order m 
with the term in square brackets divided by tt. The sec- 
ond term in brackets is proportional to the penetration 
depth in the mirror. 
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Figure 2: Schematic lateral view of a micro-cavity with opti- 
mized mirrors. The system is composed by an input-output 
waveguide, two engineered mirrors and a spacer. Each engi- 
neered mirror has coupling periods facing the waveguide and 
the spacer. 
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Figure 3: Q-factor of micro-cavities with no (circles), one (tri- 
angles) and two (squares) coupling periods in the engineered 
mirrors as a function of the resonance order. Points refer to 
numerical results, the dashed lines are the classical FP pre- 
diction and the continuous lines are the predictions from the 
FP model with effective reflectivity coefficient for the funda- 
mental mode. 



It is possible to apply this formula to our system, when 
the effective refractive index is taken equal to the calcu- 
lated effective refractive index of the fundamental mode 
of the micro-cavity. The dashed lines in Figure © are 
the Q-factors predicted by eq. Q for the micro-cavities 
with no, one and two coupling periods in the DBRs, re- 
spectively. Only the order of magnitude is got by this 
model, while the oscillations in the Q-factors found in 
the simulations are not reproduced. It is thus needed to 
go beyond this simple model and use more accurate mod- 
els which predict a non linear dependence of the Q-factor 
on the mode order. In particular the theory in Ref. 
considers that not only the fundamental Bloch mode of 
the micro-cavity is excited by the fundamental waveguide 
mode but also radiative cavity modes. What it is worth 
of mentioning is that these radiative cavity modes can 
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be coupled back into the transmitted waveguide mode, 
a phenomenon which can be described as a recycling of 
photons in the waveguide by the micro-cavity. Hence, 
not only the fundamental cavity mode contributes to de- 
termine the Q-factor of the micro-cavity but also the ra- 
diative ones. The constructive (destructive) interference 
among the cavity modes, which eventually increases (de- 
creases) the Q-factor, depends on the cavity parameters 
(physical length, modal reflectivity). To simplify the cal- 
culations, the radiative cavity modes are described by a 
single leaky cavity mode with two complex parameters, 
a complex coefficient for the coupling of the cavity leaky 
mode with the fundamental mode and a complex effective 
index of the leaky mode itself. Thus, the whole process 
can be described by an effective reflectivity coefficient for 
the waveguide fundamental mode, which depends on the 
resonance order r e fj(m). The Q-factor is then given by a 
FP-like expression with the dressed reflectivity coefficient 
r e g(m) 



Q 



\r e s(m)\ 



\r c s(m)\< 



(4) 



The continuous lines in Figure show the predictions of 
this model when the coupling and effective refractive co- 
efficients arc fit to the simulation data. The agreement in 
this case is pretty good, suggesting that the recycling of 
the leaky mode could play an important role in determin- 
ing the performance of such micro-cavities. Even if the 
fit parameters have no strict physical meaning, we have 
checked that they take values that are compatible with 
leaky mode representation: the effective index of leaky 
mode is lower than that of the cladding, the real part of 
the coupling coefficient satisfies the energy conservation. 
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Figure 4: a) top view SEM image of the cavity with 
one segment engineered mirror. b) relative transmit- 
tance of transverse electric (TE) polarized light of the 
cavity shown in a). The points shows the experimen- 
tal data and the continuous line represent the numeri- 
cal simulation. The nominal/simulation parameters are 
A =435nm/450nm and f.f .=0.23/0.244 for the coupling pe- 
riod, A =490nm/510nm and f.f.=0. 204/0. 215 for the periodic 
mirror, and 190nm/172nm for the cavity defect. 



IV. PRELIMINARY EXPERIMENTAL RESULTS 

To validate the positive role of the coupling periods 
to increase the micro-cavity Q-factors, we present some 
preliminary experimental results for a system realized by 
FIB technique, starting from Si 3 N 4 single mode waveg- 
uide. As a side-result we show that FIB technique can be 
effectively used to form ID photonic crystals with good 
optical quality without the need to go through the lithog- 
raphy and etching steps used in other processing. 

Slab multilayer waveguides were fabricated by low- 
pressure chemical vapor deposition (LPCVD) on a 2.5 
/im thick Si0 2 and it consists of the following sequence 
of Si 3 N 4 and Si0 2 : lOOnm Si 3 N 4 , 50nm Si0 2 , 200nm 
Si 3 N 4 , 50nm Si0 2 and lOOnm Si 3 N 4 [l]|. This results in 
a total core layer thickness of about 500 nm. The core was 
capped with a 500-nm-thick cladding Si0 2 layer. Lithog- 
raphy and etching defined channel waveguide geometries, 
whose nominal widths ranged from 1 to 10 /im. The 1D- 
PhC structures were defined on the channel by using a 
30 KeV Ga+ Focused Ion Beam (FIB) [HQ. In our ex- 
periment, we used the LEO-ZEISS 1540XB CrossBeam®, 
comprising a high resolution FIB column to mill directly 



in combination with a high precision scanning electron 
microscope (SEM) for precise positioning and inspection 
of the fabricated nanostructures in real time. Using an 
ion current of about 100 pA and controlling the FIB by 
a pattern generator (RAITH ELPHY) , the pattern is di- 
rectly written on the sample surface; using a total ion 
dose of 400 mA/cm 2 , a depth of about 1.3 /zm for each 
structure was obtained. Figure (0J a shows a top view 
scanning electron microscopy (SEM) image of a Fabry- 
Pcrot samples with one coupling period. 

The optical properties of integrated optical microcav- 
ities have been characterized by coupling in light from 
a tunable laser (1300-1630nm, 2mW) through a single 
mode polarization maintaining tapered fiber, mounted 
on a nano-positioning system. Two linear polarizers and 
a half-wave plate are used to control the polarization of 
the input signal. The collection system was provided by 
a near field microscope objective matched to a variable 
zoom mounted on a high-performance InGaAs infrared 
camera controlled by the LBA-500 Spiricon beam ana- 
lyzer software. A prism beam splitter allows to direct 
the transmitted signal to a calibrated photodiode (Ge 
detector) to perform intensity measurements. Normal- 
ization of the transmission intensity is done with respect 
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to a nearby reference waveguide without the ID photonic 
crystals. 

Figure Jlp) shows the relative transmittance of trans- 
verse electric (TE) polarized light for the engineered FP 
microcavity of Figure Qi). The measured spectrum 
shows a well defined band-edge at 1380 nm and the res- 
onant peak at 1515 nm. The first FP oscillation (cor- 
responding to the interference of Bloch-modes of the ID 
system) can be observed between 1300 and 1380 nm. The 
measured full width at half maximum (FWHM) of the 
resonance peak is 15nm which corresponds to a quality 
factor Q = 105. 

The spectrum is compared to the numerical calcula- 
tion (line in Fig. 0Jd) and a reasonable good agreement 
is found for the spectral position of the FP stop-band, 
the resonance position and width, while the baseline in 
the photonic band-gap region is not reproduced due to 
diffraction losses and light leakages not taken into ac- 
count in the calculations. The parameters (length of 
air slits and waveguide segments) used to reproduce the 
spectrum are slightly different from the nominal one by 
about 5%-10%. The Q-factor calculated with classical 
FP model is larger by a factor two than the measured 
one. This can be a good indication that anti-recycling 



of the leaky modes into the fundamental mode occurs in 
this sample. Clearly, other effects can affect the measured 
low Q values which are mostly related to the quality of 
the ID photonic structures, however photon recycling ef- 
fects are surely playing a role and only a more systematic 
study of these engineered structures could weigh the rel- 
ative importance of the photon recycling effect. Such a 
study is under way. 



V. CONCLUSIONS 

In this work we have presented preliminary testing of 
the effect of photon recycling which occurs in a ID pho- 
tonic crystals when guided and radiative modes interfere. 
Simulations and first measurements show that this effect 
could not be underestimated in optimizing the Q-factors 
of micro-cavities based on ID photonic crystals. 
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